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ABSTRACT 

We report on Chandra X-ray Observatory ( Chandra ) High Energy Transmission Grating spectra of the dipping 
low-mass X-ray binary 1A 1744—361 during its 2008 July outburst. We find that its persistent emission is well 
modeled by a blackbody ( kT ~ 1.0 keV) plus power law (T ~ 1.7) with an absorption edge. In the residuals of 
the combined spectrum, we find a significant absorption line at 6.961 ± 0.002 keV, consistent with the Fexxvi 
(hydrogen-like Fe) 2- 1 transition. We place an upper limit on the velocity of a redshifted flow of v <221 km s~ 1 . We 
find an equivalent width for the line of 27 + _ 2 eV, from which we determine a column density of (7 ± 1) x 10 17 cm -2 
via a curve-of-growth analysis. Using XSTAR simulations, we place a lower limit on the ionization parameter of 
>10 3 6 erg cm s _I . We discuss what implications the feature has on the system and its geometry. We also present 
Rossi X-ray Timing Explorer data accumulated during this latest outburst and, via an updated color-color diagram, 
clearly show that 1A 1744—361 is an “atoll” source. 
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1. INTRODUCTION 

1A 1744—361 is a neutron star low-mass X-ray binary 
(LMXB) discovered by the Ariel V satellite (Davison et al. 
1976; Carpenter et al. 1977). This source is a transient LMXB, 
and several outbursts have been observed with a number 
of missions, most notably the Rossi X-ray Timing Explorer 
(RXTE). In 2001, Emelyanov et al. (2001) discovered a burst 
from the direction of 1A 1744—361, however, it could not 
be unambiguously identified as a thermonuclear burst. The 
first thermonuclear (Type I) burst from 1A 1744—361 was 
discovered by Bhattacharyya et al. (2006a) using RXTE. This 
Type I X-ray burst exhibited a 530 Hz burst oscillation which 
provided a measurement of the neutron star’s spin frequency. 
The burst also provided an upper limit of d < 9 kpc, under the 
assumption that its maximum luminosity could not exceed the 
Eddington luminosity for a 1.4 M e neutron star (Bhattacharyya 
et al. 2006a). 

This source also shows regular incidents of intensity “dips” in 
its X-ray emission. According to Frank et al. (1987), such dips 
are produced by obscuring material associated with a structured 
accretion disk. Dips would only be visible for inclination angles 
60° < i < 80° (Frank et al. 1987). In principle, dips could 
occur every orbital cycle, however, because of variations in the 
obscuring material in the disk, the presence of dips often varies 
from cycle to cycle and from source to source. The dips observed 
in 1A 1744— 361, as is the case for other dipping LMXBs, exhibit 
complex structure (see Bhattacharyya et al. 2006a, Figure 4). By 
measuring the spacing between two closely spaced dip episodes 
observed by RXTE, Bhattacharyya et al. (2006a) estimated the 
orbital period of 1A 1744—361 to be 97 ± 22 minutes. 

Narrow absorption features have been seen in the per- 
sistent emission of many “dipping” LMXBs. For example, 
4U 1916—05 (Balman 2009; Iaria et al. 2006; Juett & 


3 Also at CRESST, University of Maryland Baltimore County, Baltimore, 
MD 21250, USA. 


Chakrabarty 2006; Boirin et al. 2004), XB 1624-490 (Xiang 
et al. 2009; Balman 2009; Iaria et al. 2007), 4U 1323-62 
(Balucinska-Church et al. 2009; Balman 2009; Church et al. 
2005; Boirin et al. 2005), 4U 1746-371 (Balman 2009), and 
EXO 0748—676 (Jiminez-Garate et al. 2003; Church et al. 
1998). By far, the most prominent features are the Fexxv 
(He-like) and Fe xxvi (H-like) lines. These features allow one 
to probe the structure, dynamics, and evolution of the mate- 
rial surrounding the neutron star. Thus, the spectral proper- 
ties of 1A 1744—361 are of great interest. Using RXTE data, 
Bhattacharyya et al. (2006b) found that the persistent spectrum 
of this source is well modeled by a Comptonized blackbody 
model (Bhattacharyya et al. 2006b), and they also found evi- 
dence of a broad (~0.6 keV) iron emission feature at ~6 keV 
and an iron absorption edge at ~8 keV. RXTE is not sensitive to 
narrow spectral features given its coarse spectral resolution and 
large FOV, thus, during the 2008 July outburst of this source, 
we triggered Chandra X-ray Observatory ( Chandra ) target-of- 
opportunity (ToO) observations for this purpose. 

In Section 2, we present our analysis of the Chandra 
(Section 2.1) and RXTE (Section 2.2) data from 1A 1744—361. 
In Section 2. 1 . 1 , we present high-resolution spectra of the source 
in which we find anFe xxvi absorption feature. In Section 2.2.1, 
we present an updated color-color diagram, and in Section 3 we 
discuss our results. 

2. ANALYSIS 
2.1. Chandra Observations 

During 1A 1744— 361’s last outburst in 2008 July, we trig- 
gered three Chandra ToO Observations. These observations are 
summarized in Table 1 . In order to search for narrow absorp- 
tion features such as those seen in other dipping LMXBs, the 
data were taken with the High Energy Transmission Grating 
(HETG; Canizares et al. 2005) aboard Chandra. The HETG 
consists of two arms, the High Energy Grating (HEG) and 
the Medium Energy Grating (MEG). The HEG is sensitive to 
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Table 1 

1A 1744-361 Observation Log and Spectral Properties 


Parameter 3 

2008 Jut 5 

2008 Jul 6 

2008 Jul 7 

Total 

Obs. ID 

9884 

9885 

9042 


Exposure (ks) 

24.2 

21.0 

30.4 

75.6 

Count rate (counts s -1 ) 

4.9 ± 0.02 

4.86 ± 0.02 

5.02 ± 0.02 

4.927 ± 0.009 

Continuum model 

Nu b (xlO 22 cm” 2 ) 

0 40+0-03 
u -* u -0.03 

0 35 + 0.03 
-0.03 

0 /1/1+Q-Q3 
u '**— 0.03 

0 414+0-006 
U ‘* 1 *— 0.015 

£ Ed ge c (keV) 

7 56 +0 - 08 
/ - JD — 0.05 

7 77+0.01 
/•/j-o.oi 

7 699 +0013 
' - u **-0.008 

7 645 +0 006 
/•°* -0.015 

^Edge^ 

i o+0.3 
0.2 

3 5 +1-3 

d.J-07 

3 5 +l3 
■ 3 - j - 0.6 

1 1 ' 7 + 0.10 

i - i 7 -0.09 

kT bh e (keV) 

0 no+0.08 

U-^o_0.08 

1 05 +011 
AUJ -0.08 

1 no +0-07 
1 -0.06 

1 O£+ 0-02 

1 .UO _ 0 03 

Rbb f (km) 

4 g+ 0-8 
*’°-0.7 

4 4+0-7 
*■*— 0.6 

4 4+0-4 
*■*-0.3 

4 3 +0 - 2 
*- J -o .2 

Lbol 8 (xlO 36 erg s -1 ) 

O £+0.3 
2 -6-0.3 

3 0 +0 - 6 

d-U-04 

7 c+0.5 
D-D— 0.4 

3 l +0 - 2 
J - 1 - 0.2 

r PL h 

| 64+0.06 
1,0 —0.06 

i ci +0.07 
i ' Ji — 0.05 

1 80+0-07 
EoU— o.06 

1 z:o+ 0.02 

EOo_ 0.02 

Power-law flux 1 (x 10 -9 erg cm -2 s _1 ) 

1 o£+0.04 
1 . 20-0 04 

1 7Q+0.04 
1 - 3y — 0.06 

1 14 +O.O 5 
a - 1 *-0.06 

1 ')c+0.03 
i ' ZJ - 0.02 

Total fiuxJ (x 10 -9 erg cm -2 s _1 ) 

1.399 

1.368 

1.241 

1.367 

XfioF PoF] 

0.70 [1490] 

0.69 [1333] 

0.77 [1883] 

0.83 [3407] 

Fe xxvi line properties 

Ei k (keV) 

6.97 ± 0.03 


6.97 ± 0.23 

6.961 ± 0.002 

A E l (eV) 

<50 


<229 

15.5 ± 2.4 

W E m (eV) 

33 -15 


<40 

21% 

FWHM n (kms- 1 ) 

<23192 


<5065 

1613 ± 21 

N° (xlO 17 cm- 2 ) 

l°-6 


<13 

7 4 +11 
' -*-i.o 

Notes. 


a All errors represent 1 a uncertainties. 
b Hydrogen column density. 
c Absorption edge energy threshold. 
d Maximum depth of absorption edge at energy threshold. 
e Blackbody temperature. 

I Blackbody radius as determined from the normalization of the blackbody component and assuming a distance of d = 9 kpc 
(the upper limit on the distance to the source). 

8 Bolometric luminosity as determined from the normalization of the blackbody component and assuming a distance of 
d = 9 kpc (the upper limit on the distance to the source). 

II Index of power-law component. 

1 Unabsorbed 2-10 keV flux of the power-law component. 

1 Unabsorbed 2-10 keV flux of the combined model. 
k Line energy, see Equation (1). 

1 Line width, see Equation ( 1 ). 
m Equivalent width of the line. 
n I'WHM of the line in units of velocity. 

0 Column density of the ion as determined by a curve-of-growth analysis, see Section 2.1.2 for details. 


photons in the 0.8-10 keV band with a resolving power of 
E/AE ~ 1000 at 1 keV, and the MEG is sensitive to photons in 
the 0. 4-5.0 keV band with a resolving power of 660 at 0.8 keV. 
The gratings were operated with ACIS, which enables sepa- 
ration of the individual diffraction orders, m = ±1, ±2, ±3, 
and a non-dispersed zeroth-order image is also obtained. We 
inspected all dispersed orders; however, for our analysis, we 
concentrated on the \m\ = 1 data because its effective area is 
more than an order of magnitude greater than the other orders 
for both the MEG and the HEG. Data in the different orders can 
be combined, however, they have different spectral and spatial 
responses and thus we avoided doing so, lest we smear out any 
narrow features. 

2.1.1. High-resolution Spectroscopy 


file for each arm and order following the standard Chandra anal- 
ysis threads. We then grouped the spectra so that there were no 
less than 60 counts per bin after background subtraction. We 
used XSPEC 4 v 12. 5.1 for subsequent spectral modeling. We fit 
the MEG data to a photoelectrically absorbed blackbody and 
power-law model. The MEG data was well fitted by this model, 
and there were no significant deviations. We then analyzed the 
HEG data and fitted the same spectral model with the addition 
of an absorption edge at ~7.6 keV, because of the HEG’s higher 
response to photons above 5 keV. The addition of the absorption 
edge was not physically motivated, rather it was included be- 
cause it provided a better empirical model of the continuum. For 
the HEG, we ignored bins above 8 keV because of their dimin- 
ished statistics. A photoelectrically absorbed blackbody+power- 
law+edge model fits well, see Table 1, however, the residuals 
showed a clear deviation from the model at 6.96 keV in the 
2008 July 5 (Figure 1, panels (1A) and (IB)) and 2008 July 7 


For each HETG arm and dispersed order, we extracted a 
source and background spectrum starting from the level 2 event 
file. We then created a response matrix and ancillary response 
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Figure 1. HEG order m = — 1 spectra of 1 A 1744—361. Panels (1 A), (2A), and 
(3 A) correspond to the 2008 July 5, July 6, and July 7 observations, respectively. 
The red curve represents the best-fit spectral model, see Table 1 for details. Panels 
(IB), (2B), and (3B) display the residuals, in terms of a from the model, after 
fitting for the best-fit spectral model. Notice the highly significant deviations 
from the model at 6.961 keV in panels (IB) and (3B). 

(A color version of this figure is available in the online journal.) 


(Figure 1, panels (2A) and (2B)) observations. There were no 
comparable deviations in the 2008 July 6 observation, however, 
a hint of the feature could be seen at 6.96 keV. We noticed that 
the feature was primarily in the m = — 1 and not the m = 1 
spectra, in principle the two should have comparable spectra, 
since they have comparable responses and effective areas in that 
energy range, however, dithering could affect the spectrum and 
smear out the feature, which is what we suspect happened to 
the order of m = 1 spectrum. We do not believe that the feature 
is an instrumental artifact because then it would be seen in all 
ni — — 1 HETG observations of 1A 1744—361. 

In order to increase the signal-to-noise ratio of the feature, 
we combined all three observations and their ancillary response 
files. The combined spectrum is displayed in Figure 2, and the 
6.96 keV feature can be clearly seen (Figure 2, inset). In the com- 
bined spectrum, the line is clearly resolved. Assuming that the 
line is broadened because of the thermal Doppler effects, then 
the line profile close to the peak should resemble a Gaussian: 


<Pg(E, kT) 


y/TjtAE 


exp 


(E - Ei) 


2-1 


Table 2 

Theoretical Properties 2 of the Fe xxv (He-like) and Fe xxvi (H-like) Lines 


Ion 

E 0 b 

r 

A,/ (x 10 14 s -1 ) 


(keV) 



Fexxv 

6.70010 

0.775 

5.033 

Fe xxvi 

6.96614 

0.408 

2.863 


Notes. 

a All parameters were provided by the XSTAR online line finding list 
http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html. 
b Energy of transition. 
c Oscillator strength. 
d Einstein coefficient. 


where E, is the line energy of the ion, and the line width, AE, is 
given by 

AE = EiikT/rmc 2 ) 112 , (2) 

where m ,■ is the ion mass and kT is the kinetic temperature 
(see Rybicki & Lightman 1979). Fitting Equation (1) times 
a normalization to the line profile, we obtain a line energy 
of Ei = 6.961 ± 0.002 keV and a standard deviation of 
AE = 15.5 ± 2.3 eV. The theoretical properties of the Fexxv 
(He-like) and the Fe xxvi (H-like) lines, the most prevalent lines 
in dipping LMXB spectra, are listed in Table 2. These values 
were obtained from the online line finding list 5 provided by the 
XSTAR emission modeling package (Kallman & Bautista 2001). 
The peak energy is consistent with the Fexxvi n = 2-1 tran- 
sition, and according to Equation (2), our measured line width 
corresponds to a kinetic temperature of kT = 254 ± 2 keV. 
Using the difference between the observed, £), and rest energy, 
Eq, of the line (see Table 2), we can determine the velocity of 
the absorbing material via (£)■ — Eo)/Ej = v/c, for which we 
find v < 221 km s -1 . We can only place an upper limit of the 
flow velocity because the wavelength difference we measure 
(AX = 0.001 A) is less than the absolute wavelength accuracy 
(AX = 0.006 A) of the HEG, see the Chandra X-ray Center’s 
Proposer’s Observatory Guide. 6 Given that our particular choice 
of grouping could affect the observed central energy of the ob- 
served Fexxvi line, and thereby bias our flow velocity mea- 
surements, we verified that there were no discrepancies in the 
observed energy of the Fe xxvi line between the grouped and 
ungrouped spectra. Fits to the ungrouped spectra for the two 
observations in which the line was observed, and for the com- 
bined spectrum, returned line energies that were consistent with 
those found from the grouped spectrum. Furthermore, the cor- 
responding wavelength shift of the line in the ungrouped spectra 
was all below the absolute wavelength accuracy of the detector. 


2.1.2. Curve-of-growth Analysis 

A curve-of-growth (COG) analysis allows one to relate the 
equivalent width. We, of a spectral line to the column density, 
N, of the absorbing ion. XSPEC provides a function which 
computes the equivalent width of a line, for which we obtain 
We = 27^3 eV. To verify our XSPEC-measured equivalent 
width, we also computed the equivalent width analytically, by 
assuming that in the spectral region of interest the continuum 
is dominated by the power-law component. Our analytically 
determined equivalent width, 29.8 eV, was consistent with the 
value found using XSPEC. 


2AE 2 


( 1 ) 


3 


5 

6 
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Figure 2. Sum of the three HEG m = — 1 spectra of 1A 1744—361 displayed in Figure 1. The red curve is the best-fit model (photoelectrically absorbed blackbody 
plus power law times an edge, see the text for details). Panel (IB): the residuals, in terms of o from the model, after subtracting the best-fit model displayed in panel 
(IB). Notice the highly significant deviations from the model at 6.961 keV corresponding to the Fexxvi transition. Inset: an expanded view of the residuals around 
the Fexxvi feature. The red curve is the best-fit Gaussian line profile, which has a centroid at Eq = 6.961 ± 0.002 keV with a width of A E = 15.5 ± 2.3 eV, see the 
text for details. The Fe xxv is not statistically significant, however, we fit for it in order to place an upper limit on its column density. 

(A color version of this figure is available in the online journal.) 


In order to calculate the column density, we computed 
theoretical COGs following Kotani et al. (2000) and Rybicki 
& Lightman (1979). The equivalent width, We, of the line is 
related to the optical depth, r via 

poo 

W E = / (1 - e~ T )dE , (3) 

Jo 

here r = Nscp, where N is the column density, 4> is the line 
profile shape, and 

s — itfe 1 /m i c, (4) 

where/is the oscillator strength (see Table 2), and m e and e are 
the electron mass and charge, respectively. Note that there is also 
a temperature-dependent term in Equation (4) which we ignored 
because Kotani et al. (2000) showed that the contribution of 
this term is negligible. The overall line profile, cp(E,kT ) is 
found by convolving the line profile due to thermal broadening 
(<pc. Equation (1)) and the line profile due to collisional 
broadening (7/>/), which has a Lorentzian profile, i.e., 

hy 

4tt 2 (E — Eq) 2 + (hy /2) 2 ’ (5) 

where y is the sum of the Einstein coefficients, A u j, over all lower 
energy states (see Table 2). Using the kinematic temperature 
found above, we then numerically integrated Equation (3) to 
obtain our COG, which is displayed in Figure 3. Using the error 
on the dynamic temperature and the equivalent width, we can 
obtain a confidence region for the Fe xxvi column density. Note 
that our feature appears to lie beyond the “linear regime” of the 
COG. We find V Fe xxvi = 7.4+' | l 0 x 10 17 cm" 2 for Fexxvi. 



Figure 3. COG for the Fexxvi absorption line observed in the persistent 
emission of 1A 1744— 361 : equivalent width vs. column density for a kinematic 
temperature of kT = 254 keV. The error on the COG due to the error on the 
temperature would be smaller than the thickness of the line. The horizontal 
dashed line indicates the position of our measured equivalent width, and 
the dotted lines are the upper and lower bounds on the equivalent width. 
The dash-dotted line represents the extrapolation of the “linear regime” of 
the COG. 

2.1.3. Ionization Parameter 

Although the presence of the Fe xxvi feature is clear, there 
was no discernible Fe xxv feature. However, an upper limit on 
the column density of the Fe xxv can constrain the ionization 
parameter q . 

In order to place an upper limit on the equivalent width of the 
Fe xxv transition, we fitted for such a line using our combined 
spectrum, see Figure 2 (inset). Then, using XSPEC, we calculated 
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Figure 4. Ionization parameter, £, as a function of the relative abundance of 
Fexxv to Fexxvi as determined by our XSTAR simulations, see Section 2.1.3 
for details. The vertical dotted lines represent our upper limit on the relative 
abundance of Fe xxv to Fe xxvi which places a lower limit on the ionization 
parameter of § > 10 3 ' 6 erg cm s _1 . 

the equivalent width in the same manner as for Fe xxvi. We then 
generated a COG for Fe xxv using the corresponding oscillator 
strength and Einstein coefficient (see Table 2). From our COG 
analysis, we place an upper limit on the column density of 
Fexxv of A^Fexxv < 6.3 x 10 16 cm -2 . To calculate f , we ran 
simulations of an ionized absorber using XSTAR. We assumed 
a spherical distribution of material and we used the standard 
spectral shape of a power law with a spectral index of two. We 
created a grid of ionization parameters and relative abundances 
using different initial choices for the electron density and 
luminosity. As an upper limit on the choice of luminosity, we 
used the bolometric luminosity as measured from the blackbody 
component, L < 3.1 x 10 36 erg s 1 (see Table 1). Our plot of 
| as a function of the ratio of the column density of Fe xxv to 
Fe xxvi is shown in Figure 4. Given that the relative abundance 
of Fe xxv to Fe xxvi is less than 0.09, this places a lower limit 
on the ionization parameter of £ > 10 3 6 erg cm s -1 . 

2.2. RXTE Observations 

The high time resolution data presented here were obtained 
from the Proportional Counter Array (PCA) aboard RXTE. The 
PCA is made up of five independent proportional counter units 
(PCUs). Each PCU is made up of three xenon/methane layers 
and an uppermost propane veto layer. The PCUs are sensitive to 
photons in the ~2-60 keV band binned into 256 channels. The 
PCA can time tag a photon to an accuracy of ~ 1 /xs. RXTE pro- 
vides two standard data modes. Standard-1 mode data provide 
light curves with 0.125 s temporal bins and are summed over 
all spectral channels. Standard-2 data provide light curves 
with 16 s time bins and are grouped into 129 spectral channels. 
As well as the standard modes, there are several user-selected 
modes available. The ones used here include GoodXenon and 
Event modes. GoodXenon mode provides the full temporal and 
spectral resolution. Event mode data provide high time resolu- 
tion but often at slightly lower time resolution than GoodXenon, 
and the spectral channels are grouped using different grouping 
schemes, but often preserving 64 spectral channels. 

2.2.1. Color-Color Diagram 

We created a color-color diagram using the Standard-2 
data. Each PCU has different energy response, and PCUs 0 and 


1 have lost their veto propane layer. PCU 2 is the most reliable, 
thus, in our analysis we used only photons from the top Xenon 
layer of PCU 2. For count rates below 20 counts s -1 , we used 
2048 s long intervals, for count rates between 20 and 40 we 
took 512 s long intervals, and for the highest count rates we 
used 256 s long time intervals. For each spectrum, we generated 
a model background using the FTOOLpcabackest. We defined 
our colors as follows: our soft color is the ratio of counts in 
the 2.30-3.54 keV band to that in the 3.54-5.19 keV band, 
and our hard color is defined as the ratio of the counts in the 
5.19-8.52 keV band to that in the 8.52-17.74 keV band. Our 
color-color diagram is displayed in Figure 5. In Bhattacharyya 
et al. (2006b), the authors reported that the source exhibited 
atoll-like behavior, however, a Z-source track could not be 
unambiguously ruled out. The additional RXTE data presented 
here clearly display atoll source behavior. 

2.2.2. Timing Analysis 

Bhattacharyya et al. (2006a) reported a tentative (2.3cr) 
detection of an 800 FIz quasi-periodic oscillation (QPO) for 
1A 1744—361. Using all the available RXTE event data of 
1A 1744—361, we searched for low-frequency and kHz QPOs. 
In our search, we broke up our light curves into 250 s intervals. 
We found no highly significant CS>3a) QPOs. The most signif- 
icant QPO we found was a low-frequency QPO at 29 Hz QPO 
for an observation on 2008 July 17 (Obs. ID 93155-01-06-03) 
at a significance of 3.0a. 

In the RXTE archive, we find a total of three thermonuclear 
bursts from 1A 1744—361. These bursts are displayed in 
Figure 6. The first burst occurred on 2005 July 16 (Obs. 
ID 91050-05-01-00) and was first reported by Bhattacharyya 
et al. (2006a). This burst is the brightest observed from this 
source thus far (see Figure 6, panel (a)). Using its luminosity, 
Bhattacharyya et al. (2006a) placed an upper limit on the 
distance to the source of d <9 kpc. The second and third bursts 
occurred on 2008 June 24 and 2008 July 25, respectively (see 
Figure 6, panels (b) and (c)). Bhattacharyya et al. (2006a) had 
found a significant burst oscillation in the first burst. We searched 
the other bursts for burst oscillations but found no comparable 
signals. 

3. DISCUSSION 

In this paper, we presented Chandra HETG observations of 
the dipping LMXB 1A 1744—361 during its latest outburst. 
Although the source’s persistent spectrum was well modeled 
by a simple continuum model, there was significant evidence 
for a narrow absorption feature at 6.96 keV. The feature is 
consistent with the wavelength of the Fe xxvi (hydrogen-like Fe) 
n = 2-1 transition to the wavelength accuracy of the detector. 
We find no evidence for statistically significant lines due to other 
transitions. Here, we discuss the properties of the continuum and 
line emission and we discuss what implications both have on the 
geometry of the source. We also discuss the limitations of our 
data, compare this source to other X-ray dippers, and speculate 
on what we may learn from future observations. 

3.1. The Continuum Spectrum of 1 A 1744—361 

The continuum spectra of dipping LMXBs are generally well 
modeled by two components, one due to emission near the 
surface of the neutron star, and a Comptonized component 
possibly due to an extended Accretion Disk Corona (ADC: 
see Church & Balucinska-Church 2004). Often, a blackbody 
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Figure 5. Color-color diagram of 1A 1744—361 as observed by RXTE. Left: hard color vs. soft color, see Section 2.2.1 for details. Right: hard color vs. intensity. 
From both panels notice that 1A 1744—361 exhibits atoll-like behavior. 



Figure 6. Thermonuclear X-ray bursts observed from 1A 1744—361 with RXTE. 
The light curves are background subtracted 0.125 s light curves generated using 
the full PCA energy band. The bursts were observed on 2005 July 16 (a), 2008 
June 24 (b), and 2008 July 25 (c). 


is used to characterize the central emission and a power law 
to model the Comptonized component. Such a two-component 
model provided an adequate fit for 1A 1744—361, however, 
we note that the properties of the Fexxvi were insensitive to 
our specific choice of continuum model. Using our continuum 
model, we measured a total X-ray luminosity of Lx = (1.94 ± 
0.02) x 10 36 erg s -1 in the 1-30 keV band, assuming a distance 
of cl = 9 kpc as determined by Bhattacharyya et al. (2006a). 


Church & Balucinska-Church (2004), using data from several 
dipping LMXBs, determined that the total 1-30 keV X-ray 
luminosity is directly related to the size of the ADC: 

0.88±0.16 

cm. (6) 

Thus, using Equation (6), we find that the size of 1A 1744— 36 1 ’s 
ADC is tadc = ( 1.1 ± 0.5) x 10 9 cm. According to Church & 
Balucinska-Church (2004), the radius of the accretion disk, /-disk, 
is related to the radius of the ADC by 

rdlsk = rADC (2^)’ (7) 

where P is the orbital period of the source and A t is the 
dip ingress time. Using two closely spaced dip episodes, 
Bhattacharyya et al. (2006a) were able to estimate the orbital 
period of 1A 1744—361 as 97 ± 22 minutes. 7 We found no 
additional dip episodes in the RXTE data, and only a single 
short (~2 s) dip in the Chandra data on 2008 July 5. All the 
dips had durations in the 2-25 s range. Taking 2 s as an upper 
limit on the ingress time, together with our lower bounds on both 
the orbital period and padc, we find, using Equation (7), that the 
accretion disk has a radius of rdi s k > 4.3 x 10 11 cm. Thus, for 
1A 1744—361 the ADC covers <0.3% of the accretion disk. 

3.2. The Fe xxvi Absorption Feature and Constraints 
on the Geometry of the System 

Although each individual spectrum of 1A 1744—361 was 
well fitted by a simple continuum model, the residuals showed 
clear evidence for an absorption feature at 6.961 keV in the 
first and third Chandra observations. This feature is well 
resolved in the combined spectrum. We place an upper limit 
on the velocity of a redshifted flow of v < 221 km s” 1 . We 
measure a line width of A E = 15.5 ± 2.3 eV. Under the 
assumption that the line is broadened due to thermal effects, 
the measured line width corresponds to a dynamic temperature 
kT — 254 ± 2 keV. Using this temperature, we performed a 
COG analysis for Fexxvi and determined a column density 
for Fexxvi of A^exxvi = 7 x 10 17 cm~ 2 . We do not find 
statistical evidence for a feature at the Fe xxv energy. However, 
our upper limit on the equivalent width of an underlying Fe xxv 


f ADC = 10 S 


10 35 e 


7 Note, however, that these dip episodes were separated by a data gap, thus 
not precluding a shorter period. 
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feature allows us to place a lower limit on the ionization 
parameter of £ > 10 3,6 erg cm s _1 , as determined by our XSTAR 
simulations. Furthermore, our limit on £ allows us to constrain 
the geometry of the system. According to Tarter et al. (1969), 


where L is the luminosity of the source (see Table 1), r is the 
distance from the central object to the absorber, and n e is the 
electron number density given by 


Ne_ 

Ar 


(9) 


Here, N e is the electron column density and Ar is the thickness 
of the absorbing slab. From Equation (8), we find a distance 
to the absorber of r ~ 2 x 10 11 cm, if we assume a scale 
height of A r/r ~ 1 . This distance is safely within the accretion 
disk radius found earlier. Given that this distance is only an 
upper limit, an important question is whether it is possible that 
the absorption line is emanating from within the ADC. Taking 
r = r AD c, we find Ar = 0.006r = 6.4 x 10 6 cm. For such a short 
scale height, and given the plausible inclinations for the system, 
it is hard to reconcile observing both the feature and dipping. 
Thus, it would appear that the feature is emanating outside of the 
ADC. Velocity broadening also places the absorber at >2r A DC- 

3.2.1. Cylindrical Geometry? 

Narrow absorption features have been seen almost exclusively 
in dipping LMXBs. This, together with the lack of evidence of 
orbital variability in the features, except during dipping, has led 
authors such as Boirin et al. (2004) to suggest that the absorbing 
material is distributed in a cylindrical geometry. In such a 
configuration, all LMXBs would exhibit absorption features, 
but would be observed preferentially in dippers because of their 
nearly edge-on geometry (Diaz Trigo et al. 2006; Boirin et al. 
2005). Here, we discuss what such a geometry would imply for 
1A 1744—361. If the incident angle is just grazing the top of this 
cylindrical distribution of absorbing material, then that could 
explain why the dips are so sparse and why even a small excess 
in continuum emission can render the absorption feature so 
difficult to detect. In this scenario, the height of this cylindrical 
distribution is given by Ar = r sin 9, where 60° < 6 < 80°. 
This, together with Equations (8) and (9), yields r < 2xl0 n cm. 
Note that given that we only have a lower limit on £ , the absorber 
can actually be much closer. The fact that this value for r is safely 
within the accretion disk radius (see Equation (7)), together with 
the fact that this scenario can explain, albeit in a qualitative way, 
the properties of the source, is very intriguing and worthy of 
further study with future observations. 

3.2.2. Variability in the Fe xxvi Feature? 

It is curious that the Fe xxvi feature is clearly present in all 
the Chandra observations except in the 2006 July 6 observation. 
This discrepancy cannot be due to dipping as only a single 
2 s long dip was present in all the Chandra observations — too 
short to significantly affect the total spectrum. Furthermore, the 
Chandra observation containing the dip showed evidence for 
the Fe xxvi feature. 

The properties of the continuum during the 2006 July 6 
observation may provide insight into the cause of the variability 
of the Fexxvi. The spectrum of the 2006 July 6 observation 


has a harder power-law index and a higher power-law flux 
than the other observations (see Table 1). A harder spectrum 
with a weaker Fexxvi feature, indicative of a much more 
ionized absorber, is reminiscent of the behavior exhibited by the 
microquasar GRS 1915 + 105. Neilsen & Lee (2009) found that 
the equivalent width of the Fexxvi feature in GRS 1915 + 105 
decreases as the hardness of the continuum increases. They 
suggested that the absorption features in GRS 191 5 ±105 
emanate from an accretion disk wind and that in the harder 
states the corona photoionizes the wind, making it transparent, 
thereby reducing the equivalent width of the lines (see Neilsen 
& Lee 2009). Although this interpretation is intriguing, line 
spectra due to accretion disk winds, although common from 
black hole systems, are virtually never seen in neutron star 
systems. However, a blueshifted line due to an accretion disk 
wind was observed by Miller et al. (2011) in the transient 
pulsar IGR J 17480— 2466. Also, a similar scenario has been 
suggested to explain the properties of the neutron star systems 
IGR J00291 +5934 (Paizis et al. 2005) and 1A 0535 + 262 
(Reynolds & Miller 2010). Although in Sections 3.2 and 3.2.1 
we show that even for different geometries the feature is most 
likely emanating far from the central source, unfortunately, 
we only have an upper limit on the flow velocity of the line. 
However, if we assume that the lines are emanating from a 
disk wind, we can use the line to place constraints on the 
properties of such a wind. The luminosity of the wind is given 
by L w = (l/2)/«wV 2 , where v is the flow velocity and rh \ y is the 
mass flow rate of the wind given by m w = Anr 2 n e m p v(Cl/An), 
where m p is the proton mass and Cl/ An is the covering fraction. 
We can use Equation (8) to write the wind luminosity and wind 
mass flow rate as 


L w = 2nm p v 3 (y) (10) 

and 

thw = Anm p v (11) 

respectively. Using Equations (10) and (11), we find that 
1A 1744—361 has a ratio of wind to accretion luminosity of 
Lw/Tacc < 2.8 x 10~ 5 (£2/4 jt), and a ratio of wind mass flow 
rate to mass accretion rate of m^/m. dcc < 11(Q/47t). Thus, 
even for small covering fractions, the mass flow rate of the wind 
can be comparable to the accretion rate. Interestingly, these 
fractions are in fact higher than what Lee et al. (2002) found for 
GRS 1915+ 105, 8 specifically, L w /L dcc < 7.4 x 10~ 7 (£2/47r) 
and mw/m acc < 1.3(£2/47r). With our upper limit on the flow 
velocity ( v ) and lower limit on the ionization parameter (£), 
we cannot rule out the possibility that the line emanated from 
a wind. However, notice that the ratios L w / L acc and m\v/hi. dcc 
are functions of v and / . Therefore, it is not clear why source 
shows similar behavior to the microquasar GRS 1915+105, 
while other dipping LMXBs, having comparable values for 
v and £, do not. Furthermore, IGR J 17480— 2446, which shows 
clear evidence of a wind, has L\y/L ac c ^ 2.7 x 10~ 2 (Q/47r) and 
rhw/m dCC » 52(Q/47r) (see Miller et al. 2011) — much higher 
than 1A 1744-361. 


8 Note that the accretion luminosity for GRS 1915 + 105, as a black hole 
system, is given by L acc = r]m acc c 2 , where Lee et al. (2002) assumed an 
efficient of tj ~ 0.1. For the neutron star systems 
1 A 1744—361 and 1 A 0535 + 262, we assumed a disk luminosity of 
0.5GMm acc /R for the accretion luminosity. 
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Figure 7. Error bars are the residuals after subtracting the best-fit continuum 
model from the spectrum of the 2006 July 6 Chandra observation. The smooth 
curve is the best-fit Gaussian emission line. The spectrum of the 2006 July 6 
showed no significant deviations in its residuals after subtracting the best-fit 
continuum model. However, fitting this spectrum with a continuum and line 
models line properties that are consistent with those of the other observations. 
Thus, we cannot rule out the possibility that the line was present in this 
observation as well, but rendered indiscernible because of changes in its 
continuum emission. 


As intriguing as a disk wind interpretation is, we caution that it 
is possible that this additional continuum emission at the higher 
end of the spectrum could have made any underlying Fe xxvi 
feature seem less significant. To verify this, we refit the 2006 
July 6 spectrum, including an Fe xxvi absorption feature, and 
holding the line centroid and width fixed at the values found 
from the total spectrum. The normalization of the line and 
all the parameters of the continuum were free to vary during 
the fit. The best-fit model is shown in Figure 7. The value of 
the resulting normalization of the line is consistent, within the 
lcr uncertainties, with that of the total spectrum. Furthermore, 
despite the change in the continuum, the equivalent width 
of the feature is not significantly different (within the 1 a 
uncertainties). Therefore, we conclude that we do not have the 
statistics to concretely judge whether the properties of the line 
truly changed during the 2006 July 6 observation, or whether the 
line remained the same but was indiscernible in this observation 
because of the excess continuum emission. 

3.3. Variability and Dipping 

Competing models have been put forth to explain the vari- 
ations in the properties of the continuum and absorption lines 
during dipping and non-dipping intervals. It has been proposed 
that these variations are due to progressive and partial covering 
of an extended ADC by an extended absorber (see Balucinska- 
Church et al. 2011; Church et al. 2005, 1998, for example). 
Authors such as Boirin et al. (2005) and Diaz Trigo et al. 
(2006), on the other hand, propose that an ionized absorber can 
explain dipping and non-dipping intervals, without the need of 
an extended ADC. Unfortunately, with only a single short (2 s) 
dip in our Chandra observations, we cannot compare the spectral 
properties of the continuum and the line during dipping and non- 
dipping intervals. However, we did explore whether a partially 
covered power-law model would fit the continuum as well as a 
blackbody+power law. Such a model did indeed characterize the 
data just as well as a blackbody + power law. We note that the 
properties of the line, however, were insensitive to which con- 
tinuum model we chose. Similar to the blackbody+power-law 


model, for the partially covered power law the only parameters 
that significantly changed during the 2006 July 6 observation 
were the index and normalization of the power-law component. 
Note that this observation showed no dips; therefore, if the 
properties of the feature indeed did vary during these obser- 
vations, then that is definitely unusual, as Boirin et al. (2004) 
find that the absorption features show no variability with orbital 
phase — except during dipping intervals. 

3.4. 1A 1744—361 and Other Short Period Dippers 

The plausible orbital period of 97 minutes (1.6 hr) for 
1A 1744—361 makes it the second shortest orbital period for 
a dipping LMXB that shows narrow absorption lines. The 
shortest period X-ray dipper, XB 1916—053, has a period that 
is nearly half of that of 1A 1744—361. Compared to other short 
period dippers, such as XB 1916—053 ( P — 0.83 hr) and XB 
1323-619 (P = 2.94 hr), we find that 1A 1744-361 is a 
factor of 10 more luminous than these sources (see Balman 
2009, Table 1). It may not be surprising then that we do not 
find statistical evidence for an Fexxv feature in 1A 1744—361. 
Combining Equations (6)-(8), we find that for sources with 
similar periods, a more luminous source will have a much more 
ionized absorber. Indeed, the lower limit on £ for 1 A 1744—361 
is consistent with the upper limits found for these sources (see 
Juett & Chakrabarty 2006; Balman 2009, for example). 

3.5. Implications for the 6.4 keV Emission 
Feature Observed with RXTE 

In addition to the narrow absorption feature reported here, 
Bhattacharyya et al. (2006b) discovered a broad (0.6 keV) Fe 
emission line at ~6 keV in the PCA data of 1A 1744—361. 
The fact that we now also see a narrow absorption feature at 
6.9 keV in the source’s spectrum could explain why the centroid 
of the broad emission feature was at ~6 keV rather than at 
the expected value of 6.4 keV. Bhattacharyya et al. (2006b) 
suggested that this broad emission line was actually a 6.4 keV 
line, which had flux removed by a narrow absorption feature at 
6.9 keV, and appeared at 6 keV after being binned with RXTE’ s 
course spectral resolution. Unfortunately, no broad emission 
features were observed in the individual or combined Chandra 
observations, however, this could be due to the smaller area 
of Chandra as compared to RXTE. By combining additional 
observations, one may be able to observe both the Fe emission 
feature and the Fexxvi feature in 1A 1744—361. Indeed, such 
a measurement was made for 4U 1323—62 by Boirin et al. 
(2005), using XMM-Newton. If future high spectral resolution 
observations of 1A 1744—361 reveal such an emission feature, 
and it is relativistically broadened, then such measurements can 
potentially constrain the compactness of the star. 

4. CONCLUSIONS 

Using three pointed Chandra HETG observations of 
1A 1744—361 during its latest outburst, we have discovered 
a significant 6.961 keV absorption feature in its persistent spec- 
trum. This is the first time such a narrow absorption feature has 
been observed in the spectrum of 1A 1744—361, and we identify 
this feature as an Fe xxvi line. The width of the line was found to 
be A E = 15.5 ± 2.4 eV, which corresponds to a dynamic tem- 
perature of kT — 254 ± 2 keV. We place an upper limit on 
the velocity of a redshifted flow of v < 221 km s - 1 . The line 
had an equivalent width of We = 27* 2 eV, for which we find 
a column density of A^Fexxvi = 7 ± 1 x 10 17 cm -2 via a COG 
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analysis. We place a lower limit on the ionization parameter 
of £ > 10 3 ' 6 erg cm s 1 . Using the source’s luminosity and 
the properties of the dips, we determined that 1A 1744—361 
has an ADC radius of ladc = (1.1 ± 0.5) x 10 9 cm and an 
accretion disk radius of rdi s k > 4.3 x 10 11 cm. The ionization 
parameter returned consistent results for the location of the ab- 
sorber, placing it within the accretion disk (close to the edge) 
but most likely outside of the ADC. We also argue that a cylin- 
drical distribution for the geometry of the absorbing material 
is consistent with our results. Using additional RXTE data, we 
generated an updated color-color diagram for the source, which 
confirms the suggestion of Bhattacharyya et al. (2006b) that 
this source is an “atoll” source. We report two bursts in addi- 
tion to the one reported by Bhattacharyya et al. (2006a). We 
find no evidence for a broad emission line in the Chandra data 
like the one discovered by Bhattacharyya et al. (2006b) using 
PCA data. Observations with X-ray-imaging telescopes such as 
Chandra during the source’s next outburst may reveal further 
narrow absorption features, as well as broad emission features. 
Such observations may constrain the structure and geometry of 
the system and are definitely warranted. 
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